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Abstract -The Ising model on annealed complex networks with degree distribution decaying 
algebraically as p{K) ^ K~^ has a second-order phase transition at finite temperature if A > 3. 
In the absence of space dimensionality, A controls the transition strength; mean-field theory applies 
for A > 5 but critical exponents are A-dependent if A < 5. Here we show that, as for regular lattices, 
the celebrated Lee-Yang circle theorem is obeyed for the former case. However, unlike on regular 
lattices where it is independent of dimensionality, the circle theorem fails on complex networks 
when A < 5. We discuss the importance of this result for both theory and experiments on phase 
transitions and critical phenomena. We also investigate the finite-size scaling of Lee-Yang zeros 
in both regimes as well as the multiplicative logarithmic corrections which occur at A = 5. 


Due to an abundance of natural and man-made 
network-like structures [1], phase transitions on complex 
networks are found in systems ranging from the physical 
to the sociological and form a fast developing theme of re¬ 
search nig. Investigations have centered on moments of 
the partition function and how they scale close to critical 
points. Inspired by the fundamental theorem of algebra, 
the Lee-Yang approach, on the other hand, is based upon 
the zeros of the partition function, which it interrogates di¬ 
rectly instead of through its moments [3]. We investigate 
the Ising model on complex networks through Lee-Yang 
zeros, to reveal an unexpected feature of importance for 
theoretical and experimental studies. 

For scale-free networks, the probability that a node has 
degree K decreases for large K as 

p{K)r^K-\ ( 1 ) 

Models on such networks have qualitatively similar crit¬ 
ical behavior to those on lattice systems. In particular, 
in the absence of a space dimension d, the exponent A 
determines collective behaviour. There are lower and up¬ 
per critical values of A for networks, analogous to criti¬ 
cal dimensions for lattices. At the upper critical values 
scaling behaviour is modified by logarithmic corrections, 


while above them, critical exponents assume mean-field 
values [5]. (This analogy has limitations; e.g., in complex 
networks the mean number of second neighbors decreases 
with increasing A, while for lattices it increases with di¬ 
mensionality.) Between the critical values, critical expo¬ 
nents, amplitude ratios and other universal quantities are 
A- or d-dependent. For the Ising model, as most models 
with discrete symmetry group, the lower critical values are 
Ale = 3, die = 1, while the upper ones are Auc = 5 and 
due — 4. 

From its inception, the Lee-Yang approach profoundly 
influenced modern-day statistical mechanics [6] . The idea 
is to generalize the magnetic field to the complex plane and 
consider the partition function as a polynomial in a related 
parameter [I]. Fisher extended this idea to the complex 
temperature plane [7]. The free energy is analytic in re¬ 
gions free from zeros. This is important because phase 
transitions are manifestations of nonanalyticities. These 
appear only when the density of zeros accumulate in the 
thermodynamic limit. Moments such as the susceptibility 
are directly expressible in terms of partition function ze¬ 
ros which contain all of the thermodynamic information 
accessible through standard approaches [8]. In particular, 
critical exponents are determinable through scaling of ze- 
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ros near the critical point and amplitude ratios from the 
angle at which they impact onto the real parameter axis. 

The full locus of partition function zeros contains ad¬ 
ditional information. The Lee-Yang theorem states that, 
under very general conditions, all zeros are purely imagi¬ 
nary in the complex field plane or on the unit circle after 
a change of variable. This circle theorem was proven for 
the Ising model on lattices independent of dimensionality 
d, size, and boundary conditions of the lattice and range 
of the interactions. It has since been extended to a wider 
class of regular-lattice models. (For a review see [8].) Re¬ 
cently, the Lee-Yang zeros of an Ising system were exper¬ 
imentally accessed for the first time using an approach 
which relied crucially on the validity of the circle theorem 

i- 

The validity or otherwise of the circle theorem for lat¬ 
tice systems is dimensional independent. Given the paral¬ 
lel roles played by A and d outlined above, one may expect 
an analogous statement to hold for annealed complex net¬ 
works. However, here we show that, while the theorem 
holds for the Ising model on a complex network provided 
A > Auc, it fails when A < Auc- This represents a striking 
difference between the A > Auc and A < Auc fat-tailed cases 
and an unexpected difference between models on networks 
and lattices. The finding has consequences for experimen¬ 
tal realizations of Lee-Yang zeros, and therefore for further 
investigations of phase transitions, on networks. We dis¬ 
cuss these consequences after presenting our analysis. 

The Hamiltonian of the Ising model on a complex net¬ 
work is 

= + ^ 5 ,. ( 2 ) 

l^m I 

Here, S'; = ±I is a spin variable, H an external magnetic 
field, the sum spans all pairs of N nodes, and the 

adjacency matrix Jim contains information about the net¬ 
work structure in that Jim = I if the nodes are linked and 
Jim = 0 otherwise. In general, the degree Ki = Jim 
varies for different I and may be characterized by a distri¬ 
bution p{K). For annealed networks, the links fluctuate 
on the same time scale as the spin variables. Therefore 
the partition function is averaged with respect to the link 
distribution as well as the Boltzmann distribution |10j . To 
implement this, each node I is assigned a hidden variable 
ki, the distribution of which is given by ©• The prob¬ 
ability of a link between nodes I and m is then given by 
Pim = kikm/{N{k)). It is easy to check that the expected 
node degree value is then EjiF;] = '^mPim = h- The 
above choice for pim leads to the Hamiltonian ([2|) with 
separable interaction. 

Applying the Hubbard-Stratonovich transformation, 
averaging over the spins can be performed exactly, leading 
to the representation for the partition function. 


^^lncosh[{xki + H)/T]^dx , (3) 

i 

having dropped, for clarity, pre-factors insignificant for the 
purpose of this study. 

The sum over nodes I may be rewritten in terms of a 
suitable integral over k with distribution function p{k) = 
C\k~^, in which c\ is a normalizing constant. One obtains 

m 

/» + 00 2 

Zn{T,H) = J dx, (4) 

where 

d(^) = ^ >"'“>> (± S'+f) ■ (5) 

•/N 

This partition function sets the phase diagram of the 
Ising model on an annealed scale-free network [I]l2] . There 
is no phase transition for A < 3. For A > 3 a second-order 
transition occurs at finite temperature = (/c^)/(fc). In 
the region 3 < A < 5 it is governed by A-dependent ex¬ 
ponents which, in standard notation for the thermal be¬ 
haviour of the specific heat, magnetization, susceptibility 
and the field-dependency of the magnetization, are [12] 

” = '’=1^' S'"!. « = (6) 

For A > 5 the exponents attain their usual Ising mean- 
held values a = 0, /3 = 1/2, 7 = 1, and 5 = 1/3 enhanced 
by the logarithmic corrections at A = 5. Thus the global 
variable A determines the critical behaviour in a manner 
similar to the space dimensionality for lattice models. 

We are interested in the zeros of the partition function 
dH) at the critical temperature. Substitute T = Tc in (|1]) 
and make an asymptotic expansion keeping the leading 
terms in 1/A^ and H, to Hnd 

^ cosh{hx)dx , 3 < A < 5, 

cosh{hx)dx , A > 5 . 

(7) 

Here, the H- and iV-dependencies are adsorbed in a single 
variable /i, the explicit form for which differs in different 
regions of A; 

r , 3<A<5, 

'■=] , A = 5, (8) 

[ A >5. 

Note that Z{h) is independent of A when A > 5. In par¬ 
ticular, the logarithmic corrections at A = 5 reside in 
the form of h in the middle expression of ([8]). There¬ 
fore, working in terms of the reduced variable h, the 
loci of partition function zeros for A = 5 and A > 5 
are the same. For A < 5 the values of the coefficients 
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Fig. 1: Solutions of the system of equations Re Z{h) — 0 (dia¬ 
monds, red online) and Im Z{h) = 0 (triangles, including along 
the imaginary axis) for the partition function 0 when A>5 
(left) and A = 4.5 (right). The Lee-Yang zeros are given by the 
intersections of contours of different types and are denoted by 
circles (blue online). All zeros are imaginary only when A > 5. 



Fig. 2: Lee-Yang zeros in the complex e^-plane. The diamonds 
(red online) correspond to A>5 and the circles to A = 4.5. The 
plot depicts only those zeros that fit into the region shown. 


a(A) = —c\ dyy '^(Incoshy — y^/2) > 0 can be found 
numerically and are listed in m- 

Let us now consider the partition function © in com¬ 
plex magnetic field h = Ke h + i Im h. We are interested 
in its zeros hj at the critical point. The index j identifies 
the sequence of zeros relative to the origin in the upper 
half-plane with j = 1 identifying zero closest to h = 0. In¬ 
creasing iV, the zeros approach the origin and their scaling 
at criticality is governed by the exponent a M- 




(9) 

Comparing ([8]) and (O one 

finds 


f A-2 

3 < A < 5, 

(10) 

(J ■=. }. A —1 ’ 

1 3/4, 

A > 5. 


The picture is reaffirmed in Fig. [2] which demonstrates 
that the circle theorem holds in the complex e^-plane for 
the Ising model on the annealed scale-free network for A > 
5. However it does not hold for A < 5, where only a few low 
index zeros lie on the unit circle, the rest being scattered 
on the complex plane. 

Thus, while our numerical examples hint that the Lee- 
Yang theorem holds for A > 5, we have established the 
surprising result that it does not hold when A < 5. To 
gain analytical insight, we can determine an asymptotic 
form for the partition function 0 when Re h = 0 in the 
limit of large Im h. We employ the Erdelyi lemma; if 

F{y) = [ (11) 

JO 


Therefore, a becomes A-dependent in the region 3 < 
A < 5, similar to the other critical exponents 0. More¬ 
over, the validity of the scaling relation a = (36/ {2 — a) is 
verified [M]. Also, from ([8]), a logarithmic correction ap¬ 
pears in the marginal case A = 5: Hj ~ iV“^/^(ln 
Again, the power of the logarithm complies with the cor¬ 
responding scaling relation m- 

The contours in Fig. [T] depict curves in the complex 
plane along which the real and imaginary parts of the 
partition functions o vanish at T = for A > 5 and 
3 < A < 5, exemplified by A = 4.5. Note that ImZ(/i) 
vanishes when h itself is imaginary because the partition 
function is an even function of h. The intersections of the 
different contours give the locations of the Lee-Yang zeros. 
When A > 5 all zeros plotted are on the imaginary axis. 
But for A = 4.5 only the first three zeros are imaginary 
and zeros of higher order (hj for j > 3) have non-vanishing 
real parts. This type of behaviour is found for all other 
values of A in the region 3 < A < 5; there is a finite 
number, Af, such that hj is imaginary for j < Af and has 
non-vanishing real part for j > Af. Moreover, Af decreases 
with decreasing A. The numerically determined values of 
Lee-Yang zeros are listed in Tabled] for various values of 
A < Auc = 5, and for A>5. When A = 4, A/” = 2 while for 
A = 3.5, A/'=l. 


where a > 1, /? > 0 and the function f{x) together with 
all derivatives vanishes at the upper integration limit, the 
large-y asymptotic behaviour is [16] 

CO 

E k~\- 8 

a-kV ^ , (12) 


where Ofc = exp . Putting h = ir 

one can represent the partition function for 3 < A < 5 in 
the form m- 


Z(ir) = ^ Re f x^^ ^ ^dx. (13) 

A — 1 Jo 

From ini), the asymptotic expansion is 

OO 

Z(ir) ^ (A — 1) bfcr~ ^ qq ^ (14) 

k=0 


with coefficients 


bk 


/yo) p 

5fc! 




7r(fc(A — 1) + 5) 

-h 1 

5 

COS 

10 


(15) 
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Table 1: The first five zeros hj for different values of A. 



A>5 

A = 4.5 

A = 4 

A = 3.5 

J = 1 

i3.495 

f3.495 

f3.569 

i3.762 

i = 2 

ffi.784 

i6.933 

f7.823 

1.875-bi7.212 

J = 3 

i9.636 

f9.474 

2.418-bai.466 

3.659 -b f9.496 

j = 4 

H2.229 

0.589-b 112.848 

4.014-b 114.174 

5.138-bill.351 

j = 5 

H4.650 

2.297-bH6.346 

5.446-b 116.574 

6.435-bil2.983 


The asymptotics for the A > 5 partition function © 
can be evaluated by steepest descent, leading to 

3^3 /r\i' 

2 U/ ’ 

r ^ oo , A > 5. (16) 

The trigonometric function in (1161) ensures the partition 
function alters in sign as a function of r. This signals that 
it vanishes on the imaginary axis for large (inhnite) values 
of r, suggesting that Af is infinite. If the pattern of zeros 
observed in FiglHis generic in that hj has a non-vanishing 
real part only when j > Af, this indicates that all zeros 
are on the imaginary-/i axis, suggesting that the Lee-Yang 
circle theorem is obeyed in the complex fugacity plane. 

In Fig. [3] we compare numerically calculated absolute 
values of the integral \Z(ir)\ ([7|) with its asymptotic ex¬ 
pansions m and (HU) for different A. The behaviour 
for A > 5 is qualitatively different from that at A < 5: 
whereas in the former case the function oscillates even in 
the asymptotic regime (meaning that the number of zeros 
is unbounded), this is not the case for A < 5. There, after 
a finite number of oscillations the function approaches its 
asymptotic, bounding the regime of zeros. Moreover, the 
number of oscillation decreases with decreasing A. 

In summary, besides delivering finite-size scaling, we 
found an unexpected result for the locus of Lee-Yang ze¬ 
ros for the Ising model on annealed scale-free networks, 
namely, while the circle theorem is obeyed for A > Auc = 5, 
it fails for A < Auc = 5. This is surprising because it means 
that the similarity between roles played by A (for complex 
networks) and d (for regular lattices), which control uni¬ 
versal aspects of phase transitions, do not extend to the 
fundamental Lee-Yang level. 

Recently Peng et al. related imaginary magnetic fields 
associated with a bath of Ising spins to the quantum co¬ 
herence of a probe [9]. This allowed them to identify the 
experimentally obtained times at which the quantum co¬ 
herence disappears as imaginary Lee-Yang zeros. This 
demonstration of a physical manifestation of Lee-Yang ze¬ 
ros is important at a fundamental level and points to new, 
experimental ways of studying zeros in many-bodied ma¬ 
terials. However, an identification of the type made in 
[9] relies crucially on the Lee-Yang theorem in that it can 
only be made for imaginary zeros m- The results herein 
show that while these are accessible for lattices and com¬ 
plex networks for sufficiently large A, not all network zeros 
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Fig. 3: Numerically calculated absolute values of the partition 
function along the imaginary axis: \Z{ir)\ with r = Imh. The 
asymptotic expansions using Eq. d and (d in the limit 
r —>■ oo are also plotted. The symbols “-I-” indicate Z{ir) > 0 
(black and gray online correspondingly) and circles indicate 
Z{ir) < 0 (red and pink red online). Panel (a) is for A > 5 and 
Panel (b) for A = 4.5. The number of zeros is finite for A < 5. 


are accessible in this manner for A < 5. A challenge for 
experiment is to find another way to access them. 

Finally we address the question of why the theorem is 
obeyed for some but not all values of A. It was proven 
[IB] that the theorem holds for any Ising-like model with 
ferromagnetic interactions. One may therefore anticipate 
that the Hamiltonian © for a single network realization 
should have the Lee-Yang property, since the adjacency 
matrix elements Jij are non-negative. However, averaging 
over an ensemble of networks amounts to taking a sum (or 
integral) of functions and there is no guarantee for A < 5 
that the zeros of the sum preserves the property. Indeed, 
one can consider heuristically the Landau free energy for 
scale-free networks, F{M) ^ -I- ... [T9] . 

When A < 5 the second term on the right is dominant for 
small order parameter M, delivering A-dependent critical 
exponents. The partition function is also A-dependent and 
its zeros are complex. When A > 5, however, the first term 
on the right dominates. This leads to mean-field critical 
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behaviour which is independent of A. Since mean-field 
theory obeys the Lee-Yang theorem, so too must averaged 
annealed complex networks with A > 5. 
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